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i. Introduction 

Wind blowing over a water surface generates waves in the water 
by physical processes which cannot be regarded as known. It is 
true that from one point of view the problem is merely mathematical, 
if it is accepted that the air and water motions are subject to the laws 
of classical mechanics, which are in the present instance the Navier- 
Stokes equations of viscous motion (Lamb, 1932, §328) with appro¬ 
priate boundary conditions on the moving interface. Unfortunately, 
as in the simpler problem of homogeneous turbulence (Batchelor, 
1 9 S 3 )y the present resources of mathematics are insufficient for 
a solution. Progress may be possible when suitable physical simpli¬ 
fications are made. (A typical example of such simplifications arises 
in the fluid flow past thin wings where viscosity is neglected in the 
body of the fluid, while near the wing a boundary layer calculation 
is adequate.) The real problem is then how to find the proper 
simplifications, and the recent history of fluid dynamics suggests that 
close co-operation between theory and experiment will be needed 
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for its solution. Much work in the past has been theoretical hypo¬ 
thesis without experiment, or experiment without a theoretical basis. 

As soon as one tries to think of possible simplifications, questions 
like the following suggest themselves. Can viscosity or surface 
tension be neglected ? What part is played by the air boundary layer 
and the air turbulence level ? At what distance (if any) do the wave 
motions at different points cease to interact significantly? (Here 
interaction through the air motion must be included.) What part 
(if any) is played by minute roughnesses on the sea surface ? How 
does the wind select the principal wavelength generated by it? 
In answer to such questions a number of hypotheses have been 
put forward tentatively, but before these acquire any real value 
they must be compared with observations and experiments, which 
in this subject are so difficult that they have only rarely been 
made. Some hypotheses have thus remained without proof or 
disproof for many years, and through lapse of time are now taken 
more seriously than their authors would have wished. 

The literature, then, contains many papers to which a definitive 
character cannot yet be attached, and I have not found it possible 
to fit them together into a simple coherent picture. I have there¬ 
fore thought it best not to aim at a complete survey and have 
chosen rather to treat a few topics of which a knowledge is likely to 
assist future research. Among the omissions is recent Russian work 
which appears, however, to have run nearly parallel to work in the 
West. This essay is devoted to the physical aspects of the subject, 
and little space has therefore been given to wave measurements 
without a theoretical background, although their importance in 
civil engineering is acknowledged. Most of the work described here 
deals with waves on deep water. Additional complications arise 
when there are shallow-water effects. 

The following notation is used: When co-ordinates (x>y>z) are 
used, they are rectangular with the a:- axis in the direction of the 
wind, and the #-axis vertical. The suffixes a and w refer to air and 
water respectively, the density is denoted by />, the kinematic 
viscosity by v, the surface tension by T. The wind velocity is U a ; 
I 7 m in. is the least wind velocity capable of raising waves. The wave 
height is the vertical distance between crest and trough, the wave 
velocity is denoted by c , the wavelength by A. 
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2. Calculations of initial instability 

A laminar flow is postulated on which an infinitesimal perturba¬ 
tion is superposed, and the problem is to decide whether the per¬ 
turbation grows or decays. If perturbations of every wavelength 
and frequency decay, then the flow is stable and the laminar flow 
is likely to persist; otherwise transition to turbulence may be 
expected. The following calculations apply where the air flow has 
been made laminar artificially and have no immediate application to 
sea waves since air flow in the open is always turbulent. Studies of 
the flow past a flat plate have shown, moreover, that the problems 
of initial instability and of the final turbulent wave spectrum are 
distinct, and that no information about the final spectrum can be 
derived from stability calculations. Nevertheless, such calculations 
are of interest because their predictions are precise and can in 
principle be checked experimentally, and also because a closer 
analysis of the mechanisms involved in the energy transfer from the 
wind to the perturbations may suggest useful analogies. 

2.1. Kelvin -Helmholtz instability 

Physical model . Viscosity in the air and the water is neglected. 
The thickness of the boundary layer along the interface is small 
compared with all other lengths occurring in the problem (in 
particular with the wave height) and is neglected. There is then 
a steady motion, possibly unstable, in which the interface is a 
horizontal plane, the velocity U a of the air is constant, and the water 
is at rest; the problem is to determine the stability of this motion 
when the interface is slightly perturbed. The motion is controlled 
by gravity and by surface tension at the interface. The perturbation 
is so small that the equations of motion may be linearized. 

The perturbation may be resolved by Fourier’s integral theorem 
into its harmonic components, each of which may be considered 
separately. Suppose, then, that the interface is given an infinitesimal 
deformation varying as exp [iA(# — ct)\ % where k is prescribed and 
c is to be determined. It was shown by Lord Kelvin (Sir W. 
Thomson, 1871; Lamb, 1932, §268) that 

^_ Pa ^a + | Pw~ Pa _|_ _ PwPa £/2 

Pw ~^Pa l ^ Pw Pa Pw Pa {Pw "b Pa)^ / 
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When c is real, the wave does not tend to grow and the motion is 
stable for this deformation; when c is complex, the amplitude grows 
exponentially with time, and the motion is unstable. The growth 
continues until non-linear effects become important, or until the 
assumed physical model becomes inadequate in some other way. 
For any given values of U a and k it is easy to calculate the expres¬ 
sion { }; if and only if it is positive, is the motion stable. It was 
remarked by Kelvin that { } is positive for all k if U a is sufficiently 
small, in fact, if 

U a < zigiTtp-lpaHpw +Paf ( Pw-Pa )* 

= 6socm./sec. approximately. 

When U a exceeds this critical velocity there is a band of values 
of k for which the initial deformation grows exponentially, and the 
band width increases with U a . The most unstable wavelength is 
27 T Tig-l(p w -p a )- 4 =i - 7 cm., 
and the wave velocity at the critical wind velocity is 

(P*+P^P« C/ «= r2 5 x 10-3 £/„=°- 8 cm -/sec. 

Kelvin’s theory thus predicts that no waves will grow unless the 

wind speed exceeds 650 cm./sec. Helmholtz (1868; cf. Lamb, 1932, 
§232) had suggested earlier that the motion is unstable for all values 
of U a when surface tension is neglected and the wavelength of the 
disturbance is sufficiently small. In nature waves may be observed 
at wind speeds well below 600 cm./sec. In fact, the boundary layer 
will rarely be so thin that the present theory can be applied. (For 

the connexion between theKelvin-Helmholtz theory and boundary- 

layer instability see p. 223.) 

2.2. Boundary-layer instability (Lock, i 954 > Wuest, 1949) 

The following stability calculations will be better understood after 
an examination of the stability of the flow past a semi-infinite plate. 
The calculation of this is similar to Lock’s but simpler, the results 
are more easily presented and interpreted, and the theory has been 
verified experimentally. The present survey is not the place for 
a detailed description; for an excellent account of the theory and 
its experimental verification the reader is referred to Schubauer and 
Skramstad (1947). It is suggested that during a study of Locks 
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theory the following conclusions about flow past a flat plate be kept 
in mind: 

(i) A perturbation of given frequency cannot be classified as 
stable or unstable. Instead, the theory predicts zones of decay and 
zones of amplification. A perturbation starting near the leading 
edge always decays initially; as it travels downstream it may pass 
through a zone of amplification and grow sufficiently to cause 
transition from laminar to turbulent flow; still further downstream 
the theory predicts another zone of decay if transition has not taken 
place and the flow outside the boundary layer is uniform. The zones 
of decay and amplification are separated by a curve of zero 
amplification. 

(ii) It appears to be an experimental fact that (with uniform flow 
outside the boundary layer) transition to turbulence takes place 
sooner or later in practice, even at the low levels of free-stream 
turbulence which have been attained. The transition region moves 
downstream when the free-stream turbulence is reduced, but not 
indefinitely, because of external disturbances. 

(iii) The frequency spectrum of the fully developed turbulent 
flow cannot be inferred from the stability calculation. 

(iv) Agreement between theory and experiment is obtained if, 
and only if, care is taken to reproduce experimentally the assump¬ 
tions of the theory. 

Lock's calculation; physical model . A stream of air is flowing over 
deep water in a uniform gravitational field. The viscosity in the air 
and the water causes the formation of a boundary layer (possibly 
unstable) on both sides of the horizontal interface. It is assumed 
that the boundary layer is laminar and starts at a definite point on 
the interface. (This might be realized experimentally by sucking 
off at this point the boundary layer formed upwind.) There is no 
pressure gradient in the direction of the wind, outside the boundary 
layer the air is moving uniformly, and the water is at rest. The 
steady motion in the boundary layer is first calculated (Lock, 1951)* 
and the problem is to determine the stability of this motion when 
it is perturbed by an infinitesimal sinusoidal oscillation of wave¬ 
length A in the direction of the air flow. It is assumed that the wave 
height is much smaller than the boundary-layer thickness, and that 
A is small compared with the distance x from the boundary-layer 
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origin, so that the thickening of the boundary layer in one wave¬ 
length can be neglected. Surface tension at the interface is in¬ 
cluded. 

With these limitations the linearized equations of motion have 
been reduced by Lock (1954) to a pair of ordinary differential equa¬ 
tions of Orr-Sommerfeld type (Goldstein, 1938, p. 197, eqn. (8)), 
and curves of zero amplification have been calculated. For the flat 
plate a single curve in terms of dimensionless variables is sufficient, 
but in the present problem the gravitational field introduces scales 
of length and time, and there is no advantage in the use of dimension¬ 
less variables. A separate calculation is needed for each value of U ai 
and has been made by Lock for U a = 100 to 300 cm. sec., at intervals 
of 50 cm./sec. The curves of zero amplification are presented in 
terms of dimensionless variables zirxjX and Z 7 Tv\ 0 x^jU\X\ then, al¬ 
ternatively, in terms of jc and A. The latter set is shown in fig. 1, and 
may be compared with Schubauer and Skramstad, i 947 > * 3 > for 

the flat plate. A complication is that for each prescribed set of values 
of U ay x and A two values of the (complex) wave velocity can be 
found, corresponding to Lock’s ‘air’ and ‘water waves, and the 
curves of zero amplification corresponding to these are not identical. 
(Lock remarks that the ‘water’ waves travel with a velocity nearly 
equal to the velocity of free gravity waves in water, while the ‘ water’ 
motion due to the ‘air’ waves is small. The ‘water’ waves may thus 
be more prominent in an experiment.) Each region of the (#, A) 
diagram requires two labels, indicating the stability of water and 
‘air’ waves respectively, as is shown in fig. 1. The calculation is 
very laborious, and amplification factors have not yet been calcu¬ 
lated. According to the figure, amplification takes place at x = 100 cm. 
when U a = ioocm./sec., and so in a sense a wind of ioocm./sec. is 
capable of raising waves. At this point it may be well to recall the 
conclusions relating to the flow past a plate. In the absence of 
experimental information there is no justification for supposing that 
in the present problem the limitations are less narrow than before, 
and the provisional conclusions are then, with a few obvious changes, 
the same as in the last section. The theory does not apply unless 
there is a laminar boundary layer with a well-defined leading edge, 
or a velocity profile of the same general form. The theory predicts 
decay of waves near the leading edge, with some amplification of 
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Fig. i. Neutral stability curves of A against x t after Lock (1954). 
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a band of wavelengths a few metres downwind. The precise distance 
is of small importance since turbulence is expected to set in in the 
same region. Still farther downwind the frequency spectrum of the 
turbulent wave motion becomes independent of the details of the 
amplification in the laminar region. 

A check and extension of the calculation (possibly to include 
amplification factors) would be useful. If the conditions of the 
theory are then reproduced experimentally, the comparison of 
theory and experiment may throw light on the mechanism of energy 
transfer, as was suggested earlier. In unpublished work Lock has 
studied the relation between his theory and the Kelvin-Helmholtz 
theory. He finds that the Kelvin-Helmholtz wave velocity is 
obtained when first x ^00 while the parameters (v a x/X*U a )i, 
{v w xjX 2 U^ proportional to the boundary-layer thickness remain 
fixed; afterwards the boundary-layer thickness ->o. During these 
limit operations A and U a remain constant. 

Wuesfs calculation . (The physical model is the same as Lock's.) 

Unlike Lock, Wuest (1949) does not calculate the steady laminar 
boundary layer. Instead he applies the perturbation procedure to 
three hypothetical steady velocity profiles which approximate to 
the boundary-layer profile. 

(i) In the steady motion the action of viscosity is neglected, while 
in the perturbation it is included. The steady velocity is constant in 
the air, and zero in the water; the shear stress on the interface is 
infinite. The equations for the perturbation are again of Orr- 
Sommerfeld type. This is a physically inconsistent model. 

(ii) In the steady motion viscosity is now included, and the 
boundary-layer velocity profile is approximated by a number of 
straight lines. This simplifies the calculation but, as Wuest points 
out, this type of approximation has proved unsatisfactory in other 
branches of boundary-layer theory. 

(iii) In the steady motion the curvature of the velocity profile in 
the air is now included. It is found that this modification changes 
the conclusions from model (ii) quantitatively but not qualitatively. 

The results obtained in this paper by Wuest should also be the 
same as Lock’s; any differences are due to differences in the 
mathematical treatment. Wuest’s equation (2.5) for the pressure 
differs by a viscous term from Lock (1954, equation (11)); the latter 
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appears to be correct. Consequently some of Wuest’s computations 
may contain errors. The published results are not in a suitable form 
for comparison with Lock’s. 

3. Semi-empirical theories of wave generation 

3.1. Generation by normal pressures (Eckart, 19536). 

Physical model . (This applies in the open ocean when the air flow 
is fully turbulent.) A random distribution of normal pressure (wind 
gusts) acts over the surface of a circular region of inviscid fluid of 
radius Z), the storm area. The pressures are so small that the 
equations of motion can be linearized, and the distribution of gusts 
is independent of the waves already formed. No account is taken 
at this stage of any difference in pressure between the windward and 
the leeward side of the wave crest (sheltering), and the theory is 
unrealistic in this respect. All the gusts producing the waves are 
assumed to be similar, being in effect circular pressure regions of 
the same radius L and the same duration moving over the surface 
with the constant speed U a of the wind. The gusts have already 
been acting without change for sufficiently long to make the motion 
at all points under consideration, both inside and outside the storm, 
statistically stationary in time. The size and the strength of the 
gusts are assumed to be known from observations. 

There are good reasons for considering a model of this type. 
Water is a fluid of small viscosity which transmits pressures more 
readily than shear stresses. If shear stresses are negligible then the 
powerful methods of potential theory are available for the calcula¬ 
tion, and it is natural to do this calculation first even if shear stresses 
are actually not negligible. 

Mathematically this work combines the classical theory of Cauchy 
and Poisson (Lamb, 1932, §255) with the modern ideas of statistical 
correlation analysis due to G. I. Taylor (1921) and developed by 
Eckart (1953 a), and the results are expressed in statistical terms. 
The details of the calculations are complicated. For the waves out¬ 
side the storm area it is found that at a distance r>Z), from the 
centre, the root-mean-square wave height h (from crest to trough) 
varies as r~$ } as is obvious from energy considerations, and is also 
a function of direction. In the direction of maximum wave height 
it is found that h 2 = ^oP 2 D 2 /rL, where P denotes the root-mean- 
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square wind pressure (measured in centimetres of water) at the 
centre of the idealized storm. 

This relation was used to calculate P from estimated values of 
the other variables (h = 2*8 m., D = 425 km., r = 3800 km., L = 58m.) 
for a particular storm, and it was found that P would have to be 
1 *6 cm. of water, and of the same order of magnitude as \p a U\. 
Eckart considers that o-i cm. of water is a high value for P on the 
ocean surface, although in one case 0*5 cm. was estimated; and he 
concludes that the suggested mechanism is less effective in pro¬ 
ducing waves than the actual one. (It does not appear, however, 
that reliable measurements of P have ever been made.) Since energy 
travels away from the storm in the form of wave motion it is obvious 
that the wind does work on the waves, and so there must be some 
sheltering. Eckart finds, however, that the horizontal mean stress 
in his model is only one-fortieth of that postulated by Sverdrup 
and Munk. 

Eckart ascribes the difference between his theory and observations 
partly to local correlations between wave profile and wind pressure 
(i.e. sheltering as postulated by Jeffreys; see §3.2), and partly to the 
tangential wind stress on the water which sets up a surface current 
on which rotational waves are formed. In the absence of measure¬ 
ments of P there is a gap in Eckart’s argument, but his negative 
conclusion is reasonable; the simplest model (normal pressure with¬ 
out sheltering) is likely to be insufficient. 

3.2. Generation by sheltering (Jeffreys, 1924, 1925). 

Physical model (This is applied to all stages of wave growth.) It 
is well known that the air flow past a bluff body, such as a sphere, is 
unable to follow the surface of the body but separates from the body. 
Behind the body an irregular eddying wake is formed with which the 
so-called form drag is associated. By analogy it is suggested by 
Jeffreys that the air flow over the deformed water surface is unable 
to follow the surface but separates at each crest where on the 
leeward side regions of low velocity and pressure are formed. (The 
air flow is supposed to become reattached to the water surface some¬ 
where to the windward of the next crest.) There is thus a difference 
of normal pressure between the windward and the leeward side of 
each crest, and if the wind velocity exceeds the wave velocity, energy 
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is thereby transferred from the air to the water. This is the only 
mechanism of transfer considered by Jeffreys; the effect of the 
tangential frictional drag due to molecular viscosity is found by him 
to be negligible. For a quantitative estimate of this pressure effect 
the energy balance of the wave motion is considered. The actual 
wave profile z = £(x,y> t) is for this purpose idealized into a single 
regular sine wave moving with a constant velocity c in the direction 

of the wind, £(*,y, t) = a cos-^- (x-ct)> and the normal pressure is 

then a periodic function of x which may be analysed by Fourier’s 
theorem. The only term occurring in the calculation of the work 
done on the wave is the pressure component in quadrature with £ and 
therefore in phase with 3 £/dx t which Jeffreys writes as 
sPaiUa-cysydx, 

where s is a number called the sheltering coefficient . Reasons are 
given why s should be of the order of unity. The work done by the 
wind on the water is then 

cjjpf x dxdy. Sp .(U.-c)Hjj(l$d*dy 

per unit time when U a > c . (For U a < c the modification is obvious.) 
If this exceeds the dissipation of energy by molecular viscosity, the 
waves grow, otherwise they decay. (For shallow water the dissipa¬ 
tion in the bottom boundary layer must be included; see Bondi, 
1942.) The dissipation is calculated as if the water motion were 
irrotational. 

It can now be shown that the most unstable wave is two- 
dimensional, and that the wave grows if 


<U a -cy>&gei, 


if surface tension is neglected. For given U a the left-hand side has 
a maximum when c = \U m and there is a wave (with c < U a ) which 
grows provided that 

lT u a : 


TJZ \ 4 ^w Pw & 
rU n ^ 


s Pa 

and in particular the least wind capable of raising waves is 
TT -( v wPw&\^ 


where v w =0-018 cm. 2 /sec., g = 981 cm./sec. 2 , pjp w = 0-0013; and so 
n. = cm./sec. According to observations by Jeffreys, U min . 
is about 110cm./sec., and so s is about 0-3. The corresponding 
wavelength znc 2 /g = § nU^Jg = 8*1 cm., which is beyond the range 
of capillarity. The dependence of U mhlt on s is not critical, since 
a reduction of s by a factor of ioincreases U mhlt by a factor io^ = 2*15. 

It will have been observed that a number of mathematical 
assumptions had to be made in the model which cannot strictly be 
justified, the most important being the replacement of the sea 
surface by a single sine wave. It is reasonable to suppose, however, 
that these do not give rise to wrong orders of magnitude in the 
subsequent calculations. More serious is the lack of evidence for the 
assumed physical picture of the air flow, and the energy argument 
used to derive s from U m ^ is devious (cf. §6 in the present survey). 
For a check, direct measurements of pressure over the wave profile 
are needed, and these are difficult because of the motion of the 
profile. If the moving profile is replaced by a stationary rigid profile, 
however, measurements of pressure are feasible and have been 
made on at least three occasions (Stanton, Marshall and Houghton, 
1932, hereafter referred to as Stanton; Motzfeld, 1937; Thijsse, 
1951) since the publication of the work of Jeffreys. The conclusions 
are not completely in accord. Stanton and Motzfeld agree that for 
low waves the pressure difference is too small for the required 
energy transfer, but they do not agree about the details of the 
pressure distribution; Thijsse’s experiments (published in abstract 
only) support the sheltering hypothesis but relate to motion over 
a single crest placed on a plane surface. A few details of these 
experiments are now given. 

Measurements relevant to the hypothesis of sheltering . Stanton’s 
experiments (Stantonetal. 1932). (i) The first series of measurements 
was made in a wind tunnel 50 ft. long, and of cross-section 1 by 1 ft. 
The smooth solid wave profile consisted of twenty-seven sine waves, 
of height/wavelength ratio 0-206, the wavelength increasing linearly 
with distance. The length of the 1st wave was 5-1 cm., of the 27th 
wave 21-6 cm. The mean thickness of the wave profile was 5 in. The 
normal pressure was measured through small holes on the 10th and 
27th wave for several values of U at the air flow being turbulent 
throughout. 
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(ii) Because in the first series the wave height was a considerable 
fraction of the tunnel height, the second series of measurements 
was made in a tunnel of cross-section 3 by 3 ft. The waves in this 
series were of constant wavelength, the height/wavelength ratio 
being 0-4. The mean thickness of the wave profile was 1 in. Pressure 
measurements were made on two profiles, for which A = 7*6 and 
2*6 cm. respectively. On the former, measurements made at 40A 
and 80A from the entry were in good agreement (see fig. 2); on the 
latter, measurements were made at 122 A from the entry. The static 
pressure distributions were measured through small holes in the 
model. For comparison with the hypothesis of sheltering the 
pressure must be analysed into its harmonic components: If 
£=acos 277 A:/A is the profile of the corrugated surface, then the 
coefficient of sin znx/A. in the Fourier development of pjp a U* is 
znas/k. Stanton has made Fourier analyses of five sets of measure¬ 
ments; the corresponding values of s are shown in table 1 (after 
Sverdrup and Munk, 1947). Although the waves are steeper than 
most ocean waves, the sheltering coefficient is seen to be only about 
one-tenth of the value obtained by Jeffreys (one-fiftieth in the last 
set). Stanton also notes the very interesting fact that the pressure 
variation over the waves is much less than is predicted for irrota- 
tional inviscid flow. 

Table 1 


Experiment 

A (cm.) h — 

2 a (cm.) 

U a (cm./sec.) 

s 

| 

f 10*8 

2*2 

325 

0*036 

First series, section 

1 10*8 

2*2 

470 

0*047 

1 ft. square 

| 21*6 

4*4 

330 

0*068 

1 

l 21*6 

4*4 

580 

0*090 

Second series, section 

3 ft. square 

7-6 3-0 

Table 2 

1400 

0*006 

A (cm.) 

h = za (cm.) 

No. of waves 

s 

(i) 3° 

i *5 


3 

0*034 

(ii) 3° 

3*0 


3 

0*024 

(iii) IS 

i *45 


6 

0*028 

(iv) is 

2*0 


6 

0*11 


Motzfeld's experiments (1937). Mot? ield has also measured the 
air flow over a smooth rigid corrugated surface placed in a wind 
tunnel. His wind tunnel was of height 20 cm., width 60 cm. and 


P/P m U . 
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Fig. 2. Distribution of normal pressure on the surface of solid model waves, 
according to Stanton et al. (1932), second series of experiments. O 40A from 
inlet, U a = 1400 cm. sec.; + 80A from inlet, U a = 1400 cm. sec.; x 8oAfrom 
inlet, C7 a =i400 cm./sec., model reversed; □ 80A from inlet, U a = 465 cm./sec. 


Direction of wind-► 


Profile of wave on the 
same horizontal scale 
as the pressure curves 



Fig. 3. Distribution of normal pressure on the surface of solid model waves, 
according to Motzfeld (1937). • Profile (i); O profile (ii); © profile (iii); 

O profile (iv). 


length 600 cm.; the working section was of length 90 cm. Four 
different wave profiles were constructed, the first three being sine 
waves, the fourth a sharp-crested profile consisting of circular arcs 
intersecting at 120° (table 2). The static-pressure distributions on 
the profiles were measured by a shaped static tube and also (on the 
last three profiles) through small holes in the profiles. They are 
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shown in fig. 3 and should be compared with fig. 2. The differences 
are striking. Whereas Motzfeld’s distributions are nearly sym¬ 
metrical about the crests (except for (iv)), Stanton’s are unsym- 
metrical. The air flow does not separate from the profiles, except 
over the sharp-crested profile (iv). It should be noted that the flow 
over the profiles is not independent of conditions in the entry 
length. In fact, Motzfeld has shown that significant changes are 
produced over profile (ii) by roughening the entry. Nevertheless, 
the sheltering coefficients derived from Motzfeld’s results (by his 
equation (28)), and shown in the last column of table 2, are of the 
same order of magnitude as Stanton’s, except for the sharp-crested 
profile (iv). The pressure variation along the profile is less than is 
predicted for irrotational inviscid flow along a corrugated wall, 
e.g. about one-third for profile (iii). This, too, agrees with Stanton’s 
observations. What we really need to know, however, is the character 
of the air flow when it has passed over a large number of waves, and 
it seems probable that Motzfeld’s six waves are not sufficient for 
this purpose. Stanton showed that in his experiments a final state 
was attained, and this check should always be made. The final state 
may depend on the cross-section of the wind tunnel. 

Thijsse's experiments (1951). These experiments have been pub¬ 
lished in abstract only, and only a few details can be given here. 
The smooth solid profile was of a single wave, of slightly asym¬ 
metrical shape, imitating the profile of a wind wave. The slope was 
smaller than in the experiments of Stanton and of Motzfeld. 
A stream of water was run over the wave, the pressure was measured 
through holes in the profile, and it was found that the pressure 
difference between crest and trough agreed well with the prediction 
of irrotational theory. The distribution of pressure was such that 
the energy transfer by it to a moving asymmetrical wave would be 
sufficient to account for the growth of the wave. All these conclusions 
disagree with Stanton’s and Motzfeld’s, whose experiments were 
made on wave trains and are therefore perhaps more relevant to the 
problem of wave generation than Thijsse’s. The use of water instead 
of air may also have had an unexplained effect, e.g. in the level of 
turbulence. 

A summing-up on the sheltering hypothesis . The evidence of the 
three sets of experiments which have just been described on the 
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whole favours the conclusion that the pressure differences over 
a solid profile composed of a large number of waves are an order of 
magnitude smaller than the differences postulated by Jeffreys. If 
this is so, then the hypothesis of sheltering in its simplest form is 
insufficient to explain the growth of waves, and mechanisms in¬ 
volving a drag must be considered. The substitution of a solid for 
a liquid profile should, however, be remembered. New experiments 
should be made to elucidate the differences between Stanton and 
Motzfeld, and these might also include measurements of shear 
stress for which methods are now available (Preston, 1954). The 
effect of the wind-velocity profile in the wind tunnel should also be 
investigated. As for the theoretical argument, it has been seen that 
the attempt to estimate s indirectly from E/ miT1 , leads to an erroneous 
result, and this serves as a warning of the limitations of this type of 
argument. The drag on the sea surface will be discussed in §6. The 
question arises whether this is due to a viscous skin friction, or to 
flow separation and sheltering behind small roughness elements. It 
will be seen later that there is some evidence favouring the second 
alternative. 

4. Semi-empirical systems of wave forecasting; some 

physical aspects 

This survey is not concerned, with practical methods of wave 
forecasting beyond the information which they give about the 
physics of wave generation, and this is little at present. It is true 
that the methods now in use are founded both on empirical rules 
(such as the relation of Sverdrup and Munk between wave age and 
wave steepness) and on physical assumptions, but the dependence 
of the results on the latter is not critical and little can be learnt in 
this way until the agreement between forecast and observation has 
become much closer than it is now. Although the physical basis 
of these methods is therefore still uncertain the following brief 
account is included, since some of the ideas may prove valuable. 

4.1. The system of Sverdrup and Munk (1947) 

From Stanton’s experiments (see §3.2) the authors argue that the 
sheltering effect is too small to account for the observed value of 
^min. (but cf. §6) and infer the presence of another mechanism of 
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energy transfer. They remark that the tangential wind stress does 
work on the water if there is a drift of the water in the direction of 
the wind. (This drift is erroneously identified with the Stokes mass 
transport (Lamb, 1932, §250) in irrotational free waves; in fact a 
wind wave is neither irrotational nor free.) For the tangential wind 
stress the value r = 2*6x io -3 pJJ* is taken, when U a > 500 cm./sec., 
where U a is the wind at a height of 8 to 10 m. Using the Stokes 
drift velocity and a sheltering coefficient s = 0*013 the authors find 
the energy supplied to the sea; it might be assumed that the whole 
of this energy goes into the drift currents; actually the authors 
assume that the whole energy goes into the wave motion, which for 
this purpose is idealized into a single sine wave, the ‘significant 
wave’. The downward transmission of the shear stress is not dis¬ 
cussed. The energy goes partly into increasing the wave height, 
partly into increasing the wavelength, where the partition is 
governed by an empirical result introduced at this point, namely, 
that the wave steepness 2 aj A is functionally related to the ‘ wave age ’ 
c\U a and can be regarded as known when the wave age is known. 
In support of their rule the authors have plotted a considerable 
number of measurements of wave steepness against wave age; while 
these do not lie exactly on a single curve the scatter is not too large 
for reasonable accuracy in wave forecasting. The curve chosen to 
fit the observed points is of a certain form which contains adjustable 
parameters and which the authors connect with the energy con¬ 
tributions from tangential and normal stresses. The actual energy 
partition may be very different if any of the physical assumptions 
are seriously at fault as they may well be (cf. the end of §3.2). The 
actual state of the sea, however, is more intimately connected with 
the empirical relation between steepness and wave age, and the 
forecasts of wave height based on this theory have been in fair 
agreement with observations. 

4.2. Recent developments 

In the system of Sverdrup and Munk the sheltering coefficient 
was assumed constant, and the drag coefficient was assumed inde¬ 
pendent of the wave height. It is evident that these assumptions 
can be generalized by allowing additional variations. This has been 
done by Neumann (1951, 19520,6), who has also introduced a 
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sheltering coefficient for the drag, as well as a number of other 
dimensionless quantities. Like Sverdrup and Munk he uses the 
principle that all these are functions of the wave age c/U a only, 
where the approximate functional form is deduced from observa¬ 
tions. The details are complicated. For instance, the process of 
wave generation is divided into stages, based on Neumann’s own 
visual observations. The wave age is an increasing function of 
gF/U*> where F is the fetch. When the wave age reaches a critical 
valuef depending on U a , a quasi-steady state is reached in which 
the waves break and are regenerated by energy supplied from the 
wind; at the same time longer waves arise which cover a spectral 
band around an average wave travelling with the wind speed. If the 
wind continues to blow, and fetch and duration are long enough, 
longer waves up to a limiting wave age c/U a = 1*37 are generated. 
The theoretical argument, like that of Sverdrup and Munk, is based 
only on the energy balance; this is a defect. The dissipation is due 
mainly to eddy viscosity. The reader is referred to Neumann’s 
papers for details. 

Darbyshire (1952) assumes the same mechanism as Sverdrup and 
Munk, of sheltering and tangential drag, where the sheltering 
coefficient s is taken proportional to the steepness; with this assump¬ 
tion the values of s given by Jeffreys and by Sverdrup and Munk 
can be reconciled, but Stanton’s .experiments (not mentioned by 
Darbyshire) do not agree with it. A theoretical equilibrium spectrum 
for long fetches is then derived on the assumption that the wind acts 
independently on each harmonic component of the water motion; 
while this cannot be true it is difficult to suggest anything better. 
Darbyshire does not postulate a functional relation between steep¬ 
ness and wave age; instead, he uses an empirical one between 
H(T)/T and 7 1 / U a , where T is the period. 

4.3. The spectrum of sea waves 

More recently Neumann (1953), and Pierson, Neumann and 
James (1953), have developed another approach, which does not 
utilize physical assumptions. It has always been a matter of common 
observation that the sea surface does not consist of a single wave 

t The table on p. 67 of Neumann (1951) does not agree with the formula for 
the parameter p m on p. 66. 
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train but rather of a wide spectrum, although the earlier methods 
had concentrated on the prediction of a single ‘significant wave’. 
The prediction of the complete spectrum is now attempted. If a 
wave record taken near a point in the storm area is analysed into its 
harmonic components, then the spectral wave energy dB(T) in the 
period interval (T— £d!T, T-b^dT) is by definition %gp w H 2 (T)dT t 
where H(T) is the spectral wave height. The total energy per unit 
area of sea surface near the recording instrument is 

E=\gp w ^H*{T)&T 

(Neumann, 1953, equation (10); Darbyshire, 1952, p. 302), this 
definition gives the correct value for the total (kinetic + potential) 
energy of an infinitesimal irrotational wave motion. (The energy 
spectrum of a statistically stationary process can be defined 
rigorously in terms of the correlation function introduced by G. I. 
Taylor (1921); see, for example, Batchelor, 1953, chapter 3.) Since 
the equations of motion are non-linear, and since the various wave 
components interact with each other and with the wind, any 
attempt to predict the spectrum theoretically must be regarded as 
almost hopeless, and one has to fall back on experiments. 

Now consider the wave motion which is set up by a hypothetical 
constant wind U a blowing over the whole sea and starting suddenly 
at time t = 0. The development of the motion at any later time t > o 
is then statistically the same at all points. It is assumed that the 
energy spectrum of this motion approaches a limit as £->oo, which 
may be called the spectrum of a fully developed sea. (It is also 
assumed that the same spectrum is generated by a wind U a at the 
end of an infinitely long fetch.) A similar assumption is made about 
the ‘visual spectrum 1 , i.e. the frequency distribution of time 
intervals between successive crests at a fixed position after a long 
time. Neumann has made many observations of the visual spectrum 
which he identifies with the Fourier spectrum without adequate 
justification. The suggested form for the fully developed state is 

where C is a constant of dimensions (time) -1 . From observations 

C = 8*27 x io^sec. -1 ; 
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the total energy per unit area is 

E = 3-125 xio-*Uj>, 

in c.g.s. units. In the incompletely developed state the three stages 
observed by Neumann are allowed for by assuming that for short 
fetches or durations the spectrum coincides with the fully developed 
spectrum up to a certain period depending on the fetch or duration, 
while higher periods are effectively absent. For further information 
the reader is referred to the original papers. At present the approach 
is entirely empirical, no physical assumptions about the process of 
generation have been made. 

Darbyshire (1952) has measured the pressure due to waves 
generated by local winds near a recording instrument, and has 
obtained the true Fourier spectrum by a frequency analysis of the 
pressure record. This method does not require the somewhat 
dubious identification of ‘visual* and Fourier spectra, but the 
pressure record suppresses the short periods, and the measurements 
relate to waves generated in shallow water. The notion of the fully 
developed state occurs in this work also, but the curves fitted to the 
observations differ from Neumann’s in many respects, for example, 
it is found that E varies as U 

Cox and Munk (1954) have measured the glitter of the sea surface 
and have found that the wave slopes producing the glitter are con¬ 
sistent with Neumann’s spectrum but not with Darbyshire’s (where 
short waves are suppressed). The slope distribution is Gaussian, 
and the short waves form a beam of vertical angle > 130° centred on 
the wind direction; this is wider than for long waves. The root- 
mean-square slope varies as t/f, and is tan 16 0 when U a = 14 m./sec. 

Our knowledge of the spectrum of sea waves is still very 
incomplete. 

5. Experimental studies and observations 

Among recent laboratory studies of air flow over water may be 
mentioned the papers of Stanton et aL 1932; Weinblum, 1938; 
Francis, 1949, 1951; Keulegan, 1951. Some of these are concerned 
not so much with the wave motion as with the mean drag on the 
water surface and the mean surface slope due to wind action. Some 
interesting effects have been discovered. Wuest found that the fine 
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dust circulating in a wind tunnel tended to settle on the water 
surface and thereby influenced the water motion considerably. 
Keulegan found that the addition of soap to the water inhibited the 
formation of waves by wind. Ordinarily, when the water was free 
of foreign matter, a wind velocity of 3-4m./sec. caused waves, but 
when soap was introduced into the water, and the wind velocity 
was gradually increased from zero, no waves were formed even at 
i2m./sec. The length of Keulegan’s channel was about 20 m., and 
this effect was observed at all water depths from 4 to 14*5 cm. 
Commercial household detergents had the same effect as soap. 
Keulegan finds ‘no waves* at i2m./sec., but from his measured 
velocity profiles it appears that the air motion was turbulent, and 
there must have been some disturbance of the water surface invisible 
to the naked eye (see § 6). It may be noted that according to Kelvin’s 
theory a reduction in surface tension makes the surface less stable; 
the measured effect is in the opposite direction. Francis (1951) also 
thinks that the drag on the waves is controlled by the tiny wind 
ripples. More recently (1954 a) he has blown air over the surface 
of a viscous oil, and has found that at a sharply defined critical wind 
speed small ripples appeared which were unstable and grew quickly. 
This wind speed was found to be in fair agreement with U miUit as 
predicted from Kelvin’s theory; Francis believes that Kelvin’s 
mechanism is more important here than sheltering, because of the 
high viscosity of the oil. It is not clear that this theory is applicable 
here, since Kelvin postulates a frictionless fluid and very thin 
boundary layers, and these conditions are not satisfied. Also, at 
neutral stability Kelvin predicts a wave motion in the lower fluid, 
and therefore dissipation of energy when the fluid is viscous. This 
energy cannot be supplied by Kelvin’s mechanism. 

No exhaustive laboratory study appears to have been made so far 
of wave characteristics as a function of wind strength and fetch. It 
seems probable that the wind structure will play an important part. 

The use of soap and detergents in a number of studies, for example, 
by Keulegan (1951) and van Dorn (1953), raises questions of surface 
chemistry on which much literature exists (see Adam, 1941). The 
results of Keulegan indicate that the lowering of surface tension is 
not the most important effect; in fact, a surface layer may behave 
mechanically like a gas, a liquid, or a solid (Adam, p. 315). In liquid 
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behaviour, rheological effects, such as non-Newtonian surface 
viscosity, may be expected. Solid behaviour is common; the 
damping of waves by oil (Adam, p. 105) is thought to be related to 
this. These complications are perhaps better avoided; to this end 
precautions may even be needed to keep the water surface clean. 

A very remarkable series of observations in the open has been 
made by Roll (1951) in ponds among the tidal flats at Neuwerk off 
the Frisian coast. Observations were made on the initial wave 
formation at the windward side of the ponds, and on the propaga¬ 
tion of waves up to a fetch of 7° m. The depth of the water usually 
exceeded half a wavelength. The ponds were so small that a steady 
state was rapidly attained after any change in the wind velocity. 
Even quite close to the windward edge of the pond it was possible 
to observe (i) short-crested waves of length 17 cm. arranged in 
a rhombic pattern at right angles to the wind direction, (ii) capillary 
waves (of length less than 17 cm.) with long straight crests almost 
parallel to the wind. At a certain distance from the windward edge, 
diminishing with increasing wind speed, the waves of length 17 cm. 
split up suddenly into gravity waves and ripples; this is attributed 
to a transition from laminar to turbulent flow in the air boundary 
layer. The steepness of these waves may approach 1 in 7, and 
diminishes with increasing distance downwind. The wavelength of 
the gravity waves increases as they travel downwind. When the 
wind speed at a height of 35 cm. exceeds 5 m./sec. there is a change 
in the appearance of the sea surface farther downwind; below 
5 m./sec. the capillary waves are stable and the surface always 
appears smooth, while beyond that limit the surface is rough. Roll 
attributes this to a change from laminar to turbulent flow in the 
water boundary layer. The measurements of wavelengths as a 
function of wind velocity and fetch form a single consistent picture 
into which the measurements of Scott Russell (1844), Jeffreys and 
Stanton fit without effort. While not strictly relevant to this review, 
such agreement is rare in this subject, and Roll’s diagram is included 
here. 

The observations of wave velocity as a function of wavelength 



agree with the irrotational theory and also with earlier measurements 
by Wuest and Stanton, for wavelengths exceeding 10 cm. It was 
also found that the mean shear stress (measured from wind velocity 
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profiles) is in general a decreasing function of U a1 except in the 
neighbourhood of t/ a = 6m./sec. (measured at a height of 10 m., 
equivalent to about 4*5 m./sec. at 35 cm.). This velocity is also 
associated with the change in appearance of the surface (see also 
Munk, 1947). To explain the sudden transition from smooth waves 
to gravity waves and ripples Roll considers the air flow. The smooth- 
water region near the windward edge is compared to the laminar 
region near the leading edge of a flat plate, and it is found that the 
transition on the water surface takes place at a fetch F , where 
UaF/Va -3 x z ° 5 approximately, and this agrees with the usual 
Reynolds number for transition to turbulence on a flat plate. This 
model does not explain the waves observed by Roll in the smooth 
region, and here the discussion is not complete. Actually in this 
region the air flow is not comparable with the flow past a plate. For 
over the land the wind is certainly turbulent even if the roughness 
is small; on reaching smooth water the wind cannot suddenly 
become laminar, but since the roughness over smooth water is less 
than over the land a layer of reduced turbulence starts from the 
edge of the pond and thickens downwind.f Meanwhile the rough¬ 
ness of the water surface increases downwind owing to wave forma¬ 
tion, and the resulting air flow is more complex than Roll's picture. 
The wave motion in this region is probably due to the low-turbulence 
layer. Roll's achievement in difficult conditions is astonishing, but 
it seems unlikely that much more can be done in this way without 
control, or at least measurement, of the turbulent air flow. 

6 . A note on U min 

It has often been observed that strong winds raise waves more 
easily than weak ones. Also, the Kelvin-Helmholtz theory for an 
ideal fluid predicts that for small values of U a the waves do not grow 
under the action of the wind, and that a critical value (= 650 cm./sec. 
approximately) must be exceeded before they do. From these two 
pieces of evidence the idea has arisen that in real fluids also there is 
a least wind-velocity (here denoted by U^n.) capable of raising 
waves, but this is now seen to be improbable. Suppose that the 
air stream is turbulent, as will usually be the case in the open. Then 

•f* Laboratory experiments on the air flow near the join of two plates of 
different roughness have been made by Jacobs (1939). 
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random stresses are continually exerted on the interface; if the air 
motion is, for example, statistically stationary in time, a stationary 
water motion will thereby be generated. There seems to be no 
reason why this should not be described as a wave motion. The 
amplitude may be so small that the motion is invisible to the naked 
eye, but then U mhLt depends on the quality of the microscopes which 
happen to be available. The only remaining possibility is that the 
air motion is laminar. This can be dismissed, since the discussion of 
Lock’s theory and of the flow past a plate shows that this condition 
can only be brought about artificially, and that transition to tur¬ 
bulence takes place after the wind has passed over a distance 
depending on the wind velocity and the turbulence level. 

A table of measured values of U m ^ is given below: 

E /min , (cm./sec.) Source 


40 

Roll (1951) 

85 

Scott Russell (1844) 

IIO 

Jeffreys (1924) 

200-300 

Stanton et al. (1932), Weinblum, 
(1938), Wuest (1949) 

300-400 

Keulegan (1951), without soap 

790 

Francis (1949) 

> 1200 

Keulegan (1951), with soap 


The differences in this table are probably partly due to differences 
in wind structure. The theoretical arguments and the experimental 
results both point to U m ^ = o as the most probable value, provided 
that the fetch is long enough. The foregoing remarks do not exclude 
the possibility that there is a turbulent wind velocity (appropriately 
defined) at which the waves begin to grow rapidly. 


than of waves. It is known that the mean air flow in the turbulent 
friction layer over rough ground is given by an equation of the form 


U a (z) = %[s75 togio Q + 2 '5- B i(^)] > 

when the thermal stability is neutral! Here = (r 0 /p a )^, where r 0 is 
the mean tangential stress on the ground; z 0 is a length characteristic 
of the roughness of the ground; when u%z Q jv a < o* 1, the flow is said 
to be aerodynamically smooth; when u^zjv a > 3 it is aerodynamic- 
ally rough. B x is a function which is effectively zero for z Q u#/v a > 3, 
and which depends somewhat on the type of surface roughness 
(Hama, 1954, figs. 10 and 11). The experiments of Nikuradse 
(see Ellison’s paper in this volume, p. 406) relate to sand roughness; 
those of Hama to artificial random roughness similar to paint 
roughness, and to wire screens. While in the region of transi¬ 
tional flow the differences are considerable, in the fully rough regime 
the differences due to the type of roughness are fairly small. 
Hama quotes evidence that the logarithmic profile (established by 
Nikuradse for pipe flow) also applies in a boundary layer. In smooth 
flow the drag is the viscous skin friction on the plate, while in rough 
flow the drag is a form drag due to the separation of the flow behind 
individual protuberances, and is independent of viscosity. For our 
purpose it is important to note that the flow can be characterized 
as rough, transitional, or smooth from measurements in the air> not 
necessarily very near the ground, provided that the velocity profile is 
logarithmic. The velocity % is given by the slope of U a (z) against 
5’75 logio*, and then z 0 u%/v a can be found since 


575logi„^ 


UM 

»* 


5*75 l°gio 


2f 0 «* 



7. The roughness of the sea 

It has been seen that the force of the wind on the water is not 
adequately described as a normal pressure distribution acting on 
a smooth-wave profile; the tangential stress must therefore be taken 
into consideration. The generation of waves and the wave spectrum 
are governed not only by the mean drag but probably also by the 
variation of the drag between crest and trough; but so little is known 
about this that only the mean drag can be considered here, although 
its study belongs more properly to the subject of ocean currents 


where the left-hand side is known when 1% is known, and the right- 
hand side is a known function of z Q u^jv a if B x is assumed known. 
Ellison has emphasized that measurements should always be 
analysed in terms of % and z Qy and not in terms of wind velocities 
at a given height. These results apply to air flow over the ground, 
but they may be expected to hold over water, where the roughness 
may be expected to vary with wind speed, and allowance may have 
to be made for the wind drift in the water. The results of measure¬ 
ment should be expressed in terms of characteristic quantities 
16 


BDS 
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analogous to 1% and z Qi but unfortunately it is not obvious what 
these are, since the distribution of eddy viscosity in the water 
depends on the action of the wind and is not known. (If the eddy 
viscosity in the water were kzf ( T 0 lp w ), then Ellison’s calculation 
on p. 417 shows that the earth’s rotation reduces the ultimate 
thickness of the friction layer in the water to ( pJPw )* times the 
thickness of the layer in the air. The drift current thus extends 
to a depth of the order of 30 m.) When a logarithmic wind 
profile is observed, however, it is reasonable to assume in the 
absence of more precise information that the constants of the profile 
have the same physical significance as over land, in particular that 
the friction velocity 1 % = V( T o/Pa) ma Y be found from the slope, and 
that z 0 as determined from the equation at the beginning of this 
section is proportional to the size of the roughnesses causing the 
mean drag. It may be necessary to measure wind velocities relative 
to the speed of travel of the roughnesses. The question with which 
we are here concerned is the degree of aerodynamic roughness of 
the water surface. Unfortunately, the measurement of velocity is 
not an easy one, and the evidence is not yet conclusive. Thus Roll 
(1948) measured the wind in the lowest 2 m. over the shallow sea 
at Neuwerk and found that 

U a (z) = %[s-7S lo gio^ +1-85]. 

from which z 0 u#/v a = 0-63, and the flow is transitional.f It is seen 
that z 0 decreases as % increases, an unexpected result; also, z Q was 
found to be independent of the wave height. The order of magnitude 
of Zq is 3 x io -3 cm., although in the absence of waves a much smaller 
value was observed. 

Many observers think that the flow is fully rough at sufficiently 
high winds. If this is so, then viscosity is unimportant, and the 
only dimensionally correct relation is u%=gz 0 times a constant. 
Ellison (this volume, p. 409) has quoted some observations to 
support such a relation, and more evidence has been collected by 
Charnock (1955). On the other hand, an analysis of Keulegan’s 
results (1953, figs. 16 and 17) shows that the flow is rough and that 



t Nikuradse’s Bj is used here (see p. 406). Roll gives z Q u#/v a = o- 48, but 
the difference is unimportant. 
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Zq is a decreasing function of u%. Keulegan finds that z Q is of the 
order of 0*5 cm. No measurements were made to show that the 
velocity profile is logarithmic. Van Dorn (1953) has made some 
wind measurements near a pond at three different heights. Although 
these are sufficient to show that the wind profile is not exactly 
logarithmic they provide evidence that z 0 is of the order of o-1 cm., 
and that z 0 is a slowly increasing function of u*. A new set of wind- 
profile measurements would be very welcome. If it then appears 
that the flow is fully rough for high winds, this would provide 
definite evidence that for such winds the tangential drag on the 
mean wave profile is due to a sheltering effect behind fairly small 
roughnesses (of height 30^) on the water surface. This conclusion 
would presumably apply not only to the mean drag but also to the 
local drag at all points on the water surface. (Compare, for 
example, Francis, 1951, and the references in Francis, 19546.) 

Surface-slope measurements . The determination of the drag from 
wind-velocity profiles depends on the validity of the equation for 
U a (z) at the beginning of this section. There is another method 
which appears to be more direct. When wind blows over a lake, it is 
observed that wa^er piles up at the downwind end of the lake, and 
this effect has been much studied because of the damage caused by 
flooding (for references see Francis, 1954 b ). It is clear that the surface 
stress plays an important part, and if it is assumed that the wind 
stress is balanced by the component of gravity acting along the sur¬ 
face, then the inclination a of the surface is given by tana = r 0 !gp w D , 
where D is the depth, and no allowance is made for the bottom stress 

or the return current. When these are allowed for, tan a = J * 

SPu 

(Keulegan, 1951), where T b is the bottom stress ( = |r 0 for laminar 
flow, but possibly much smaller in turbulent conditions, Francis, 
i 953). Unfortunately, the effect is a small one. (Typical values of 
mean slope are 5 x1c -4 , Keulegan, 1951; io -5 , Keulegan, 1953; 
2 x io -5 , van Dorn, 1953.) By making very slight changes in the 
derivation one can obtain different equations for the slope. This will 
now be shown by a reconsideration of Keulegan’s treatment (1951), 
which is among the most careful in the literature. Keulegan assumes 
that a tangential drag acting on the surface produces a slow steady 
laminar flow, and that the total mass crossing any vertical plane is 
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zero. (This is due to the action of the end walls.) As the slope is so 
small, the pressure at any level is nearly hydrostatic: 

P — gpJP Q 4- a - z) 4- local barometric pressure. 

(The symbols are explained in fig. 5.) If the barometric pressure is 
assumed constant (any variation is easily allowed for), then 

dp _ da 


dx 


=gP< 


'dx* 


The error in this equation may be estimated from the equation of 
vertical motion; it appears that it is sufficiently accurate. The 



|d»+(r, + r») 


equation of rate of change of horizontal momentum is then 
(Keulegan, 1951, equation (15)) 

_0 

dx~ 

= -gP W D-£+(T 0 + T b ), (i) 

and if the term on the left is negligible we have the 

usual formula. Keulegan also shows, by neglecting the vertical 
motion, that / ~ \ 

u=u Hd)’ (2) 

where F(Q= 3£ 2 —2£. An important question is whether the term 
on the left of (1) is indeed negligible. Equation (2) is exact only when 
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the depth is exactly constant; if there is a small slope a better 
assumption satisfying the equation of continuity is 

U 0 D 0 ( z \ 

U ~D 0 + a'\D 0 +aJ’ 

and the vertical velocity is then 

... U 0 D 0 z d a t z \ 
W -(D 0 +afdx r \D 0 + a)’ 

from the equation of continuity. It is then easily shown that 

j: 


* 2 . _ 2 „ U 2 D 2 0 

P w u dsr i6Pw Dn + a > 


1 0 "0 

and the term on the left of equation (1) is 


~TSPv>, 


UqDq da^ 2 TJ 2 — 
V TbPw u 0 . 



(D 0 +a) 2 d» ‘ 16r '"’“ u dx 

The ratio of this to the hydrostatic term on the right is -£sUHgD 0 
approximately, where U 0 is the surface velocity. In laminar flow 

and the ratio is thus —— > an( i Keulegan 

0 6 v w dx’ 270 v% \dxj 

da o^cm. , , 

(fig. 7) shows that, when Z> 0 = 8cm., -, and the ratio 

v 6 u u cbc 1900 cm. 

is 1*2, which is not small. Thus the allowance for a small slope of 

1 in 4000 has introduced a large term which has always hitherto 

been neglected. 

Another difficulty arises from the end walls. The calculation 
shows that the horizontal momentum transport per unit time is 
approximately p w UlD^i-zp w gD% which must be balanced by 
an increased pressure on the end wall. If this were brought about 
by a local rise in level near the wall, and if the pressure were hydro¬ 
static near the wall, the rise in level would be of the order of centi¬ 
metres, and would exceed the effect of wind drag. If the mean slope 
is calculated from water levels on the shores of a lake, and not 
from measurements in the middle of the lake, then large errors 
may result in the calculated values of the wind drag. It is clear 
that the momentum transport and the pressure forces cannot be 
neglected. 

The momentum transport and the rise of level near the shore 
need not be directly associated with a surface drag. Free water 
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waves, for instance, transport momentum in the direction of pro¬ 
pagation (Lamb, 1932, §250) and cause a radiation pressure when 
they are incident on a wall, with a consequent rise in level. When 
wind blows over a lake it generates both waves and drift currents, 
and the resulting momentum is likely to be a complicated function 
of the wind strength. In many cases, e.g. when flooding of the shore 
is studied, a measurement of the rise in level may be sufficient but 
it is clearly not sufficient for the study of wave generation. 


8 . Conclusion 

The problem with which this review has been concerned is the 
physical process of wave generation by wind. Certain pieces of 
evidence have been presented which are regarded as likely to throw 
light on the problem. The present state of our knowledge is pro¬ 
foundly unsatisfactory. Nevertheless some general conclusions have 
emerged, among which the following seem specially deserving of 
notice: 

(i) Little more will be learnt by measuring water motions which 
are generated by blowing air over water in an uncontrolled manner. 
The present review has shown in several places (in particular in the 
measurements of U min ) that a precise control of air and water 
motions is needed; support for this view comes from both theory 
and experiment. Theory and experiment should be closely 
integrated. 

(ii) Energy arguments of the type introduced by Jeffreys should 
be recognized as insufficient without other confirmatory evidence. 
The experiments described in §3.2 have shown that the physical 
picture obtained in this way can be misleading. 

(iii) Some remarks in this essay have been critical, but the 
problems which emerge do not appear to be insoluble and give hope 
for the future. For instance, it should be possible to elucidate by 
new measurements the discrepancy between Stanton and Motzfeld 
and the distribution of normal pressure; or to establish whether the 
wind-velocity profile is indeed logarithmic, whether the air flow is 
aerodynamically rough, and what are appropriate parameters 
analogous to and z 0 for the combined air-water systems. This 
requires a reconsideration of the eddy viscosity in the water, and 
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leads on to a study of shear-stress transmission, and thence to 
characteristic wavelengths. Another useful approach would be 
through Lock’s problem combined with an analysis of the energy 
transfer. A study of the energy transfer near the upwind side of 
a pond would also be useful. 

My thanks are due to Dr A. A. Townsend who has freely given 
his help and advice throughout; and to Dr M. S. Longuet-Higgins 
for much helpful discussion and criticism. 
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